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Abstract 

In this paper we study the Lyapunov stability and the Hopf bifurcation 
in a system coupling an hexagonal centrifugal governor with a steam engine. 
Here are given sufficient conditions for the stability of the equilibrium state 
and of the bifurcating periodic orbit. These conditions are expressed in terms 
of the physical parameters of the system, and hold for parameters outside a 
variety of co dimension two. 
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1 Introduction 



The centrifugal governor is a device that automatically controls the speed of 
an engine. The most important one, invented by James Watt in 1788 — Watt 
governor — , is regarded as the starting landmark for the theory of automatic 
control. The historical relevance of this device as well as its importance for 
present day theoretical and technological control developments — going from 
steam to diesel, gasoline engines and electronic governors — have been widely 
discussed by MacFarlane [6], Denny [2], Fasol [3] and Wellstead - Readman 
[T2] among others. 

The centrifugal governor design received several important modifications 
as well as other types of governors were also developed. From MacFarlane 
[6], p. 251, we quote: 

"Several important advances in automatic control technology were made in 
the latter half of the 19th century. A key modification to the flyball governor 
was the introduction of a simple means of setting the desired running speed 
of the engine being controlled by balancing the centrifugal force of the flyballs 
against a spring, and using the preset spring tension to set the running speed 
of the engine" . 

This paper is devoted to the study of the dynamic stability and simplest 
bifurcations of the system coupling the hexagonal centrifugal governor with 
a spring — called Hexagonal governor — and the steam engine. See Fig. [1] 
for an illustration. The system coupling the Hexagonal governor (resp. Watt 
governor, with no spring and with vanishing horizontal edges of the hexagon) 
and the steam engine will be called simply the Hexagonal Governor System 
(HGS) (resp. Watt Governor System (WGS)). The stability analysis of the 
stationary states and of small amplitude oscillations of this system will be 
pursued here. 

The first mathematical analysis of the stability conditions in the WGS was 
due to Maxwell [7] and, in a user friendly style, likely to be better understood 
by engineers, by Vyshnegradskii [TT]. A simplified version of the WGS local 
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stability based on the work of Vyshnegradskii is presented by Pontryagin [8]. 

The oscillatory, small amplitude, behavior in the WGS has been asso- 
ciated to a periodic orbit that appears from a Hopf bifurcation. This was 
established by Hassard et al. in [4], Al-Humadi and Kazarinoff in [1J and, in 
a more general context, by the authors in [HJ HP] - 

In [10], restricting ourselves to Pontryagin's system of differential equa- 
tions for the WGS, we carried out a deeper investigation of the stability of 
the equilibrium along the critical Hopf bifurcations up to codimension 3, hap- 
pening at a unique point at which the bifurcation diagram was established. 
A conclusion derived from the properties of the bifurcation diagram implied 
the existence of parameters where the WGS has an attracting periodic orbit 
coexisting with an attracting equilibrium. 

The results of the present paper extend in a different direction the analysis 
in [9], as described below. 

In Section [2] we introduce the differential equations that model the HGS 
illustrated in Fig. [U The stability of the equilibrium point of this model 
is analyzed and a general version of the stability condition is obtained and 
presented in the terminology of Vyshnegradskii (Theorem 12.21 and Remark 
12.31) . The codimension 1 Hopf bifurcation for the HGS differential equations 
is studied in Section [3j An expression which determines the sign of the first 
Lyapunov coefficient is obtained (Theorem 13.11) . Sufficient conditions for the 
stability of the bifurcating periodic orbit are given. Two pertinent particular 
cases (no spring and vanishing horizontal edge) are calculated and illustrated. 
See Theorem 13.61 Fig. [2] and Theorem I3.8[ Fig. [3j 

Concluding comments are presented in Section HI 
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2 The Hexagonal governor system 



2.1 Hexagonal governor differential equations 

The HGS studied in this paper is shown in Fig. [TJ There, cp £ (0, f ) is the 
angle of deviation of the arms of the governor from its vertical direction axis 
Si, Q G [0, oo) is the angular velocity of the rotation of the engine flywheel 
D, 9 is the angular velocity of the rotation of Si, I is the length of the arms, 
m is the mass of each ball, H is a sleeve which supports the arms and slides 
along Si, T is a set of transmission gears and V is the valve that determines 
the supply of steam to the engine. 




Figure 1: Hexagonal centrifugal governor — steam engine system. 



The HGS differential equations can be found as follows. For simplicity, 
we neglect the mass of the sleeve and of the arms. There are four forces 
acting on the balls at all times. They are the tangential component of the 
gravity 

— mg sin<y2, 

where g is the standard acceleration of gravity; the tangential component of 
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the centrifugal force 

m(L + I sin p)9 2 cos tp, 

2L > is the distance AA' in Fig. [H the tangential component of the 
restoring force due to the spring 

— 2kl{l — cos p) sin ip, 

21 is the natural length of the spring and k > is the spring constant; and 
the force of friction 

—blip, 

b > is the friction coefficient. 

From Newton's Second Law of Motion, using the transmission function 
9 = c Q, where c > 0, one has 

2-k r^l 2r^\ • 2M + Hig . 6 . 

(p — c —U cos <£> + \- c \l ] smp cos 99 sin <p p. ( 1) 

£ \m / ml m 

The torque acting on the flywheel D is 

/ (l = 11 cos ip - F, (2) 

where / is the moment of inertia of the flywheel, F is an equivalent torque 
of the load and \i > is a proportionality constant to represent the torque 
due to the steam which decreases with the angle (p. See [H], p. 217, for more 
details. 

From Eqs. ([1]) and ([2]) the differential equations of our model are given 

by 
dp 

dip 2 L 2 , f 2k , 2r^ 2kl + mg ______ b 



—— = c — ft, cos p + h c f2 sin 99 cos <p sin p ip 

dr I \m J m ' m 

— = - {{icosp-F) (3) 
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where r is the time. 

The standard Watt governor differential equations as presented in Pon- 
tryagin [8], p. 217, are obtained from by taking L = and k = 0, 

d tp 

d7 = * 

dip 2 2 • 9 ■ b 

— — = c \l simp cosip — - snip ip (A) 

dr 'I'm 

— = - fa cos ip - F) 

Performing the following changes in the coordinates, parameters and time 
/ ml \ 1//2 , / ml \ 1//2 

x = ^y = \^n — ^ z = c \^un — n ' 

\2kl + mg J \2kl + mg J 

(2kl + mg\ l/2 L 2kl 

b ( ml \ 1 / 2 cli f ml \ _ F 
" -r ~n~i > Z 3 



m \2kl + mg J I \2kl + mg J \i 

where p > 0, < re < 1, e > 0, a > and < (3 < 1, the differential 
equations (EJ) can be written as 

x' — — = y 
dt y 

y' = ~r = P -2 2 cos x + (z 2 + re) sin a; cos x — sin x — e y (6) 

' dz ( m 

z = — = alcosx — p) 
dt y HJ 

or equivalently by 

x' = /(x,C), (7) 

where 

/(x, C) = (y, p z 2 cos x + (z 2 + re) sin x cos x — sin x — e y, a (cos x — /?)) , 
x= (x,y,z) e(0,^xlx [0,oo) 

and 

C = ((3, a, e, p, re) G (0, 1) x (0, oo) x (0, oo) x [0, oo) x [0, 1) . 
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2.2 Stability analysis at the equilibrium point 

) have only one admissible equilibrium 



The HGS differential equations 
point 



(1-^)^(1-^)1/4 

P = (x ,y ,z ) = I arccos/3,0,— 



(■ 



/3V2( p+ (1-/32)1/2)1/2 

The Jacobian matrix of / at P has the form 

/ 1 0\ 



(8) 



Df (P ) 



where 



and 



-u% -e £ 

\ -a(l -P 2 ) 1/2 I 

1 (1 _ ^2)3/2 + _ ^3) 
/3( P + (1-/32)1/2) 



(9) 



(10) 



£ = 2/3 1 / 2 (l - /3 2 ) x / 4 (l - rc/3) 1/2 (p + (1 - (3 2 )) 1/2 . 

For the sake of completeness we state the following lemma whose proof 
can be found in [8], p. 58. 

Lemma 2.1 The polynomial L(\) = poX 3 +PiX 2 +P2X + P3, Po > 0, with real 
coefficients has all roots with negative real parts if and only if the numbers 
Pi,P2,P3 are positive and the inequality p\Pi > P0P3 is satisfied. 



Theorem 2.2 If 



e > E r 



2a^ 2 (l - (3 2 fl\l - k(3) x ' 2 { P + (1 - /3 2 )i/2)3/2 



11) 



(1 - /9 2 ) 3 / 2 + P {1 - K/3 3 ) 

then the HGS differential equations |7J) have an asymptotically stable equilib 
rium point at Pq. If 

< £ < £ r 



then P n is unstable. 
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Proof. The characteristic polynomial of Df (Pq) is given by p(\), where 

-33(A) = A 3 +pi X 2 +p 2 X + P3, 

(l-/3 2 ) 3 / 2 + p(l-^ 3 ) 
Pl £ ' P2 /?(p+(l-/9 2 )V2) 

and 

2a^ 2 {l - (3 2 f/\l - Kf3) l /\p + (1 - /3 2 ) 1 / 2 ) 3 / 2 
P3_ /3(p+(l-/? 2 )V2) 

The coefficients of — p(A) are positive. Thus a necessary and sufficient con- 
dition for the asymptotic stability of the equilibrium point P , as provided 
by the condition for one real negative root and a pair of complex conjugate 
roots with negative real part, is given by ffTTT) . according to Lemma [27TT 



Remark 2.3 In terms of the HGS physical parameters, condition fill) is 
equivalent to 



where 



— V >1, (12) 
m 



/3 2 ) 3 / 2 + p - P^np 



2/? 3 / 2 (i - /? 2 ) 3 / 4 (i - nf3y/ 2 {{i - f3 2 yi 2 + P fi 2 

is the non-uniformity of the performance of the engine which quantifies the 
change in the engine speed with respect to the load (see J2J/, p. 219, for more 
details). Eq. [W\l can be written in terms of the original parameters of the 
HGS, but this expression is too long to be put in print. 

The rules formulated by Vyshnegradskii to enhance the stability of the 
system follow directly from (112j) . In particular, the interpretation of (1121) 
is that a sufficient amount of damping — b — must be present relative to 
the other physical parameters for the system to be stable at the desired 
operating speed. Condition f|T2|) is equivalent to the original condition given 
by Vyshnegradskii for the WGS (see [8], p. 219). 
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3 Hopf bifurcation analysis 



In this section we study the stability of Po under the condition 

e = e c , (14) 

that is, on the Hopf hypersurface which is complementary to the range of 
validity of Theorem 12.21 

3.1 Generalities on Hopf bifurcations 

The study outlined below is based on the approach found in the book of 
Kuznetsov [5], pp 177-181. 

Consider the differential equations 

x' = /(x,/i), (15) 

where x G M 3 and \i e W 1 is a vector of control parameters. Suppose 0151) 
has an equilibrium point x = x at fi = fiQ and represent 

F(x) = /(x,/i ) (16) 

as 

F(x) = Ax+ \ 5(x,x) + \ C(x,x,x) + 0(||x|| 4 ), (17) 
2 o 

where A = / x (0,/io) and 

j,k=i 

C • (x ' y • z) = £ 1 ^pt*l^ t ' ! ''■ 2, • (19) 

for i = 1, 2, 3. Here the variable x — xo is also denoted by x. 

Suppose (x , Ho) is an equilibrium point of (fT5]) where the Jacobian matrix 
A has a pair of purely imaginary eigenvalues A2,3 = ±itOo, u; > 0, and no 
other critical (i.e., on the imaginary axis) eigenvalues. 
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The two dimensional center manifold can be parametrized by w G M 2 = 
C, by means of x = H(w, w), which is written as 

H (w, w) = wq + wq + — — h jk w^w k + 0(|u;| 4 ), 

' j\k\ 

with hjk G C 3 , hjk = hkj- 

Substituting these expressions into f|T5l) and f|T71) one has 

H w (w, w)w' + Hw(w, w)w = F(H(w, w)). (20) 

Let p,q & C 3 be vectors such that 

3 

Ag = iu;o<?, A T p = -iu p, (p,q) = ^2piqi =1. (21) 

The complex vectors hy are to be determined so that equation ( |20l) writes 
as follows 

w' = iu w + - G 2 iw\w\ 2 + 0(|w| 4 ), (22) 

with G 2 i G C. 

Solving the linear system obtained by expanding (1201 . the coefficients of 
the quadratic terms of ( TT6l) lead to 



h 11 = -A- 1 B(q,q), (23) 

/i 20 = (2iu I 3 -A)- l B(q,q), (24) 

where J3 is the unit 3x3 matrix. 

The coefficients of the cubic terms are also uniquely calculated, except 
for the term w 2 w, whose coefficient satisfies a singular system for h2i 

(iu h - A)h 21 = C(q, q, q) + B(q, h 20 ) + 2B(q, h n ) - G n q, (25) 

which has a solution if and only if 

(p, C(q, q, q) + B(q, h 20 ) + 2B(q, h u ) - G 2l q) = 0. 
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Therefore 

G21 = (p, C(q, q, q) + B(q, (2zu; / 3 - A)~ 1 B(q, q)) - 2B(q, A^B^q, g))). (26) 

The first Lyapunov coefficient l\ is defined by 

i 1 = J_ReG 21 . (27) 
2 uq 

From ff22l its sign decides the stability, when negative, or instability, when 
positive, of the equilibrium. 

A H op f point (x ,/io) is an equilibrium point of (fl5|) where the Jacobian 
matrix A has a pair of purely imaginary eigenvalues A2,3 = Hojq, u > 0, 
and no other critical eigenvalues. At a Hopf point, a two dimensional center 
manifold is well-defined, which is invariant under the flow generated by (fl5l) 
and can be smoothly continued to nearby parameter values. 

A Hopf point is called transversal if the curves of complex eigenvalues 
cross the imaginary axis with non-zero derivative. 

In a neighborhood of a transversal Hopf point with l\ 7^ the dynamic 
behavior of the system (fT5l) . reduced to the family of parameter-dependent 
continuations of the center manifold, is orbitally topologically equivalent to 
the complex normal form 

w' = (7 + iu)w + liw\w\ 2 , (28) 

w G C, 7, uj and l\ are smooth continuations of 0, u and the first Lyapunov 
coefficient at the Hopf point [5]. When l\ < (l\ > 0) a family of stable 
(unstable) periodic orbits can be found on this family of center manifolds, 
shrinking to the equilibrium point at the Hopf point. 

3.2 Hopf bifurcation in the HGS 

From (JTj) write the Taylor expansion (fl71) of /(x). Define 

"* = (29) 
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and 



a 



1 - K/3 



(30) 



Thus, with ojq given in Eq. ffTUl) . 

/ 1 



A 



-UJr 



o \ 

a/3 1 /2 u;i 
/ 



(31) 



\ -afl 1 / 2 ^ 
and, with the notation in (fTTj) we have 

F(x) - Ax= (0,F 2 (x) + O(||x|| 4 ),F 3 (x) + O(||x|| 4 )), (32) 

where 

^ + g^ + 1 + (3 - 7 f» (1 -^ J »- 

bp 

pV 2 a(P 1/2 ui) 1/2 (P + 4/3 x/ V) x 2 ^ + (2/? 2 - 1 - p^VO xz 2 , 



and 



s (x) = -- a (3 x 2 + \ a (3 l/2 Ui 



x 



2 ' 6 
From (I3TT) the eigenvalues of A are 

Ai = -e c , A 2 = i w , A 3 = -i u . 

The eigenvectors g and p satisfying (12 1H are respectively 



^'0 



and 



P 



z cj - ie c e c u (e c + iuj ) 



2 , 2(uj 2 + e 2 c ) , 2a[3 1 / 2 Lo 1 (uj 2 + e 2 
The main result of this section can be formulated now 



(33) 

(34) 
(35) 
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Theorem 3.1 Consider the family of differential equations The first 

Lyapunov coefficient at the point for parameter values satisfying p^ ) is 
given by 

' l(fta -"- K) = - 4fe^;( £ ; + 5^ + 4^) ' (36) 

where 

R(P, a, p, k) = £ 2 c ujfal(3pa 2 - 4)(e 2 + Awl) + Pe 2 c u^l(e 2 c + 4u;g) + 
Sa/? 11 / 4 ^ 4 ^ 2 ^ 2 + 4u, 2 ) + 8a 2 /3 3 / 2 e c Vu, V - 
4a 3 /3 13 / 4 e c ^ 1 11/2 (e 2 + Wlo 2 ) + Sa 3 /? 21 / 4 ^^ 1 / 2 ^ + lOu, 2 ) + 
4aWc4(e 2 + 8a; 2 ) + 8a 3 ^e.auf 2 (e 2 + 10a, 2 ) - 
4a 3 /3 17 / 4 e cP aa; 1 15/2 (e 2 + lOw 2 ) + 4a 4 /? 6 u;fcc + 8u, 2 ) - 
2a/3 7/4 e c cja; 2 a;i /2 (14u; 4 + 3e 2 c tul(pa 2 - 1) + 30^w 2 (/)ct 2 - 1)) + 
S(i ll/2 {Aa 2 e 2 c a 2 uj 2 Q uj\ + a 4 pujl(e 2 c + 8w 2 )) - 
4a/3 15 / 4 e c au;i /2 (-14^ - 30ufal (pa 2 - 1) + 
a 2 e 2 c uf{l + p 2 ) + w 2 w 2 (3e 2 (l - pa 2 ) + 10a 2 w 2 (l + p 2 ))) + 
2a 2 ^ /2 e 2 c pujlujl(e 2 c + 4(wg + p<r 2 u;?)) + 2a0 )/A e c paujluj\ /2 
(e 2 c {ul - Zu 2 x (pa 2 - 1)) - 10u 2 (uj 2 + Zu 2 x {pa 2 - 1))) - 
2a 2 /3 4 w 2 u; 4 (2e 4 - 8u 2 (lo 2 + 3w 2 (pa 2 - 1)) + e 2 c (9uj 2 + w 2 (3 + 13,0a 2 ))) + 
f3 2 e 2 c uj 2 {e 2 c (ut - Au 2 cu 2 (pa 2 - 1) + 2a 2 u, 4 ) + 2(4u, 6 + ufa 2 (pa 2 - 1) + 

LU 2 Luf(4a 2 + 9(/)ct 2 - l) 2 ) + 8a W)) + 
2a 2 /? 7/2 w 2 c^8/xj 2 (^ + 3ujf(pa 2 - 1)) + 
e 2 (-16a 2 a; 2 + 3pVu; 2 - + 3w 2 ))). 



Proof. The proof depends on preliminary calculations presented below. 
From ([I7D, (USD, (USD and (p} one has 

S(x, y) = (0, B 2 (x, y), -a/? x x j/0 , (37) 

where 

£ 2 (x, y) = -3 1/2 u x (l - pa 2 ) x x y x + 2(3(p + u x p 1/2 ) x 3 y 3 + 
2a{^ 2 u 1 f/ 2 {{2(3 2 - I) - p^l 2 ) 
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(x x 1/3 + x 3 yi), 



C(x, y, z) = (0, C 2 (x, y, z), a/? 1/2 u;i x x y x z x ) , 

where 

C 2 (x, y, z) = 1 + (1 ~ P ^ )(3 ~ 7/32) X! m z x + 2(2/? 2 - 1 - p/^V) 
{x m z 3 + x m z 3 + x 3 y 3Zl ) - 2f3 l/2 a{(3 1/2 uj l ) 1 ' 2 {p + 4/3 1/2 ^) 

{XiViZz + Xi?/32;i + X 3 yiZx). 

Referring to the notation in (1371) . fl38|) and fl34l) one has 

B(q,q) = (0 : B 2 (q,q),aP), 

where 



(38) 



(39) 



S 2 (g,g) 



2a 2 /?(/? 1 / 2 u; 1 ) 3/2 (p + /9 1/ V) + 



^(/^vv) 3 / 2 

3(1 - pa 2 )a; 2 (/3 1/2 cu 1 ) 1/2 + iAaau Q u x {l - 2(3 2 + /^Vi) 
B(q,q) = {0,B 2 {q,q),-aP), 



(40) 



where 



B 2 (q,q) 



pl^Lolipo 2 - 1) + 2a 2 /^V(p + pWui)) 



C(9, 9, 9) = (0, C 2 (g, g, g), -ia/^V) , 



(41) 



where 



Ca(g,?,?) 



— z 



cj 2 (4 - 3pa 2 + 7/3 2 (pa 2 - 1)) + 2a 2 /? V 2 



(2/? 2 - 1 - f3 1/2 puj 1 ) - i2a/? 2 (/? 1 / 2 u; 1 ) 1 / 2 (p + 4/? 1 / 2 a; 1 ) 

The first Lyapunov coefficient is given by (|27[) . From (|35p and (|4ip one has 

-1 



Re(p,C(g,g,g)) 



2a(3 9 ^au;y i /2 (p + 4p 1 / 2 u 1 ) + 



2(3u 2 (el + u 2 ) 

+e c (^ 2 (3pa 2 - 4 + 7(3 2 (1 - pa 2 )) + /3u% + (42) 

2a 2 /? 2 a; 1 2 (l-2/3 2 + /3 1 /V 1 ) 
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From (J3S]), (JH7D, ([Ml) and 

Re(p, 2B(q,h u )) 



one has 



s c uj$u 2 (e 2 c + a; 2 



4a 4 /3 13/4 p 2 ^i 6 + 8a 3 /? 4 e c a^ 13/2 + 8a 4 /? 15 / 4 po; 7 - 4a 3 /? 5 / 2 e c pao; 1 1E ' " + 



4a 3 e c pao; 1 11/2 (2/3 7 / 2 - /3 3 / 2 ) + 

15/2 

ia^^cuf + P llA e 2 A{ul + 3o; 2 (pa 2 - 1)) + 
2a 2 /3 7/4 pc^ 3 (a; 2 + 3a; 2 (pa 2 - 1)) + 
2a 2 /3 9 / 4 u; V(a; 2 + 3o; 2 (pa 2 - 1)) - 
2a£ c aa> 2 a; 1 5/2 (2a; 2 + 3a; 2 (pa 2 - 1)) - 
2a/3 1/2 £ c paa; 2 a; 1 7/2 (2a; 2 + 3a; 2 (pa 2 - 1)) - 



(43) 



4a/3 2 e c ao; 1 5/z (a 2 a; 4 (l + p 2 ) - 3o; 2 o; 2 (pa 2 ~ 1) 



2a; 4 ) 



From ([35]), ([37D, ([34D and ([24D one has 



Re{p,B(q, h 20 )) 



where 

a, p, k 



Aai3 11 / 4 aJ 1 /2 (e 2 J 



(44) 



+ 



%a 2 pl 2 e c o 2 ujluj\ + 4a 3 apul 1/2 (e 2 c - 2u 2 )(2(3 17 ^ - /3 9 / 4 
Aa^^uf + 8a 3 /3 19 /Vu; 3/2 (£ 2 - 2o> 2 
Sa 2 ^ 2 e c a 2 pW u; 7 - 4a 3 /3 13 / 4 apo; 1 15/2 ( £ 2 - 2a; 2 ) + 
4a 4 /? 5 e c u\ + %a 2 p 9 l 2 e c ujl{a 2 pul - 4aW ) + 
2a 2 p 3 e c tu 2 uj 4 (3uj% + o> 2 (19pa 2 - 3)) 
2a 2 p 5/2 e c u 2 Luf(l6a 2 uj 2 + 3p 2 a 2 uf + 3p(o> 2 - uf) 
2a/3 3/4 aa> 2 o> 5/2 (2a; 2 (a; 2 - 3a; 2 (pa 2 - 1)) + e 2 c (Au 2 + 3a; 2 (pa 2 - 1)) 
2a{J 5/A p(TLU 2 LU 7 1 /2 (2uj 2 (uj 2 - 3a; 2 (pa 2 - 1)) + e 2 c (Au 2 + 3u 2 (pa 2 - 1)) 

-4a; 4 - 3a; 2 a> 2 (pa 2 - 1) + a 2 a; 4 (l + p 2 ) 
2a; 2 (a; 4 - 3a; 2 a; 2 (pa 2 - 1) + cAj 4 (1 + p 2 )) 
f3e c u 2 (e 2 c u J 2 u J 2 1 (u 2 + 3(pa 2 -l) 
2a; 2 (3cj 4 (pa 2 - 1) + 9cj 2 a; 2 (pa 2 - 1) + 8a 2 pa 2 a; 4 )) 



15 



Substituting (H2j|. ^3]) and into Q2ED and (|2ZD, the theorem is proved. 



Remark 3.2 The denominator of the first Lyapunov coefficient given by Eq. 
(ESP is positive. Thus the sign of h is determined by the sign of the function 
R, the numerator of l\. 

The expression for l\ depends only the parameters a, (3, p and k, although 
in the expression in $36)) appear also uj ,uji,a and e c . This is due to the 
fact that these last parameters are functions of the previous ones as shown 
in $TQ), $MB and 

Proposition 3.3 Consider the family of differential equations (Ej) regarded 
as dependent on the parameter e. The real part, 7, of the pair of complex 
eigenvalues verifies 



l'{e c ) 



2K + el) 



< 0. 



(45) 



Therefore, the transversality condition holds at the Hopf point. 



Proof. Let X(e) = A2,3(e) = 7(e) ± iio(e) be eigenvalues of A(e) such that 
7(e c ) = and to(e c ) = u , according to (1531) . Taking the inner product of p 
with the derivative of A(e)q(e) = X(e)q(e) at e = e c one has 



dA 

^' de 



q 



1 \e c )±o J \e c ). 



Thus the transversality condition is given by 

dA 



As 



y(e c ) = Re (p, 



dA 

de 



de 



q = (0, -u; ,0) , 

c 

the proposition follows from fl35|) and a simple calculation. 



(46) 
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The full expression of l\ in terms of the parameters a, j3, p, n seems too 
long to be of use in qualitative arguments. Two special cases are considered 
below for the sake of illustration. 

3.2.1 The case p = 

Corollary 3.4 Consider the case where p = 0. Then the equilibrium point 
Po in is given by 

P = (x , ?/o, Zo) = ( arccos (3, 0, Q - kJ ) , (47) 



the Hopf hypersurface (TO) is given by 



e c = e c (P, a,n) = 2a [3 3/2 (1 - k/3) 1/2 (48) 

and the numerator of l\ in (EfJ)) is given by 

a, k) = -3 + 5k(3 - (a 2 - 5)(3 2 + n(a 2 - 7)/? 3 - 

2a 2 n 2 {3 A - (a 4 - 2aV)/? 6 + a 4 K/3 7 . (49) 

If G\ is different from zero then the family of HGS differential equations (OJ) 
has a transversal Hopf point at P for e c = 2 a (3 3 / 2 (1 — k/3) 1//2 . 

Proof. The proof is immediate by substituting p = into Eqs. ([H]), f|T4"j) and 
A sufficient condition for being a Hopf point is that the first Lyapunov 
coefficient l\ ^ 0, since the transversality condition is satisfied by Proposition 
But from fT49l it is equivalent to G\ ^ 0. 



Remark 3.5 The expression shows that the "running speed" of the sys- 
tem depends monotonically decreasing on k, which is monotonically increas- 
ing on k, according to |3]]. This corroborates analytically the quotation of 
MacFarlane in the Introduction. 
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Equation (149 p gives a simple expression to determine the sign of the first 
Lyapunov coefficient fl36l) for the case p = 0. The graph Gi(j3,a, k) = 
is illustrated in Fig. [2J, where the signs of the first Lyapunov coefficient 
are also represented. The surface li = divides the hypersurface of critical 
parameters e c = 2 a /3 3 ^ 2 (I — k/3) 1 / 2 into two connected components denoted 
by S and U where h < and l\ > respectively. In Fig. [2j the (3 coordinates 
at the reference points B\ and B2 are 0.7746 and 0.5272, respectively. 




Figure 2: Signs of the first Lyapunov coefficient for p = 0. 



The following theorem summarizes the results in this subsection. 

Theorem 3.6 Consider the case where p = 0. // (f3,a,n) G S U U then 
the family of differential equations (0|) has a transversal Hopf point at Pq for 
e = e c . If (/3, a, k) G S then the Hopf point at P for e = e c is asymptotically 
stable and for each e < e c , but close to e c , there exists a stable periodic orbit 
near the unstable equilibrium point Pq. If (ft, a, k) G U then the Hopf point 
at Pq for e = e c is unstable and for each e > e c , but close to e c , there exists 
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an unstable periodic orbit near the asymptotically stable equilibrium point P$. 
See Fig\E 



3.2.2 The case k = 

Corollary 3.7 Consider the case where k = 0. Then the equilibrium point 
Po in (El) is given by 

( (1 _ 02)1/4 \ 

P = (x ,yo,zo) = (^arccos (3, 0, ^ _ ^1/2)1/2 J ' ( 50 ) 
the Hopf hypersurface [Lfy is given by 

e c = e c (P, a, p) = p+(1 _ / g2 ) 3/ 2 ( 51 ) 

and the numerator of l\ in kMfy is given by 

G 2 (f3, a, P ) = -2a 4 /3 22 + 2a 4 /3 20 (8 + 7p((l - (3 2 ) 1 ' 2 + 3p)) + 3(-2 + 
p(-9(l - /3 2 ) 1/2 + p(-15 - 10(1 - p 2 ) 1/2 p + 3(1 - (3 2 ) 1/2 p 3 + p 4 ))) - 
2/3 18 (-5 + q 2 (1 + p((l - Z? 2 ) 1 / 2 + bp)) + a 4 (28 + p(50(l - P 2 ) 1/2 + 
7p(22 + 5p((l - /3 2 ) 1/2 + p))))) + /? 16 (-86 + a 2 (16 + p(27(l - (3 2 ) 1/2 + 
2p(52 + 9(1 - /3 2 ) 1 / 2 p + 7p 2 ))) + 2a 4 (56 + p(153(l - (5 2 ) 1 ' 2 + 
7p(69 + p(33(l - (3 2 ) 1 ' 2 + p(35 + 3(1 - (3 2 ) l ' 2 p + p 2 )))))) - 
/3 14 (-41(8 + (1 - (3 2 ) 1/2 p) + 2q 2 (28 + p(64(l - fi 2 ) 1 ' 2 + p(215 + 
2p(40(l - P 2 ) 1/2 + p(34 + (1 - P 2 ) 1/2 p)))) + a 2 (70 + 
p(260(l - (5 2 ) 1/2 + p(p(630(l - f3 2 ) 1/2 + 840 + p(700 + 
p(140(l - /3 2 ) 1 / 2 + p(56 + (1 - P 2 ) l/2 p))))))))) + /3 12 (-13(56 + 
21(1 - /3 2 ) 1/2 p + 4p 2 ) + a 2 (112 + p(319(l - (3 2 ) 1/2 + p(976 + 
p(552(l - (3 2 ) 1/2 + p(494 + 51(1 - /3 2 ) 1 / 2 p + 8p 2 ))))) + 2a 4 (56 + 
p(265(l - (3 2 ) 1/2 + p(875 + p(910(l - (3 2 ) 1/2 + p(1050 + 
p(350(l - H 2 ) 1 ' 2 + p(154 + 11(1 - (3 2 ) l ' 2 p + p 2 )))))))) - 
/3 10 (-259(4 + 3(1 - f3 2 ) 1/2 p) + p 2 (-305 + 9(1 - (3 2 ) 1/2 p) + a 2 (140 + 
p(480(l - /3 2 ) 1/2 + p(1376 + p(1011(l - (3 2 ) l/2 + 2p(458 + 
93(1 - /3 2 ) 1/2 p + 23p 2 ))))) + 2a 4 (28 + p(162(l - j3 2 ) 1 ' 2 + p(546 + 

p(735(l - /3 2 ) 1 / 2 + p(875 + p(420(l - (3 2 ) 1 / 2 + p(196 + p + (52) 
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3p(9(l - /? 2 ) 1/2 ))))))))) + /5 8 (-245(4 + 5(1 - (3 2 ) 1/2 P ) + 2 (-745 + 
29(1 - f3 2 ) l,2 p + 27p 2 ) + a 2 (112 + 0(457(1 - (3 2 ) 1/2 + 0(1264 + 
0(1121(1 - P 2 ) 1/2 + 2/9(502 + ,0(143(1 - /? 2 ) 1/2 + 2/9(20 + (1 - 
P 2 ) l/2 P))))))) + 2« 4 (8 + p(55(l - (3 2 ) 1 ' 2 + 0(189 + 0(315(1 - (3 2 ) 1 ' 2 + 
p(385 + 0(245(1 - /3 2 ) 1/2 + 0(119 + 25(1 - /3 2 ) 1 / 2 /> + 3/9 2 )))))))) - 
/3 6 (-77(8 + 15(1 - (3 2 ) l/2 p) + p 2 (-970 + 3(1 - (3 2 ) l/2 p + 
112/9 2 ) + a 2 (56 + 0(272(1 - /3 2 ) 1/2 + 0(750 + 0(795(1 - /? 2 ) 1/2 + 
p(710 + 0(240(1 - P 2 ) 1/2 + p(68 + 5(1 - /3 2 ) 1/2 /9))))))) + 2a 4 (l + 
0(8(1 - (3 2 ) 1/2 + 0(28 + 0(56(1 - (3 2 ) 1 ' 2 + P (70 + 0(56(1 - (5 2 ) 1 / 2 + 
0(28 + 8(1 - (5 2 ) l ' 2 p + 2 )))))))) + /3 4 (-248 + 0(-651(l - (3 2 ) 1 ' 2 + 
0(-71O + 0(-75(l - I3 2 ) 1 / 2 + 0(143 + 27(1 - (3 2 ) 1/2 p)))) + a 2 (16 + 
0(93(1 - P 2 ) 1/2 + 0(264 + 0(345(1 - (3 2 ) l/2 + 135(1 - (3 2 ) 1/2 p 2 + 
3(1 - /3 2 ) 1 / 2 4 + 400(8 + 2 )))))) - /3 2 (-58 + P (-203(l - (3 2 ) 1 / 2 + 
0(-277 + 0(-88(l - ^ 2 ) 1/2 + 0(58 + 47(1 - (3 2 ) l/2 p + 9p 2 )))) + 
2a 2 (l + 0(7(1 - (5 2 ) 1 ' 2 + 0(21 + 0(35(1 - f3 2 ) 1 ' 2 + P (35 + 
0(21(1 - /? 2 ) 1/2 + P(7 + (1 " P 2 ) 1/2 P))))))))- 

If G*2 is different from zero then the family of HGS differential equations (Ej) 
has a transversal Hop f point at P for e = e c . 

Proof. The proof is obtained by substituting k = into Eqs. (jSJ), f|T4|) and 
( |36l) . The long expression above, being a challenge to hand calculation, has 
been performed with Computer Algebra. In the site [13] has been posted the 
main steps of the long calculations involved in this substitution. This has 
been done in the form of a notebook for MATHEMATICA 5 |14j . A sufficient 
condition for being a Hopf point is that the first Lyapunov coefficient l\ ^ 0, 
since the transversality condition is satisfied by Proposition 13.31 But from 
( |52|) it is equivalent to G 2 ^ 0. 

■ 

Equation (|52|) gives an expression to determine the sign of the first Lya- 
punov coefficient fl36|) for the case k = 0. The graph G 2 (/3,o;,0) = is 
illustrated in Fig. [31 where the signs of the first Lyapunov coefficient are also 
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represented. The surface l\ = divides the hypersurface of critical parame- 
ters s = e c into two connected components denoted by S and U where l\ < 
and l\ > respectively. At point B\ the p coordinate is 0.0478. See Fig. [31 

a: -t 




Figure 3: Signs of the first Lyapunov coefficient for At = 0. 

The following theorem summarizes the results in this subsection. 

Theorem 3.8 Consider the case k = 0. // (/3, a, p) G S U U then the family 
of differential equations (TJP has a transversal H op f point at Pq for e = e c . If 
a, p) G S then the H op f point at P for e = e c is asymptotically stable and 
for each e < e c , but close to e c , there exists a stable periodic orbit near the 
unstable equilibrium point Pq. If ((3, a, p) G U then the H op f point at P$ for 
e = e c is unstable and for each e > e c , but close to e c , there exists an unstable 
periodic orbit near the asymptotically stable equilibrium point P . See Fig\^ 
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4 Concluding comments 



In this paper the original stability analysis due to Maxwell and Vyshnegrad- 
skii of the Watt Centrifugal Governor System — WGS — has been extended 
to the Hexagonal Governor System — HGS — where a more general force, due 
to the spring, acting on the sliding sleeve of the governor has been considered. 

In Theorem 12.21 we have extended the stability results presented in Pon- 
tryagin [8] to include this more general system. See [9] for another possible 
extension. 

Concerning the bifurcations of the HGS, this paper deals with the codi- 
mension one Hopf bifurcations in the Hexagonal governor differential equa- 
tions. The general expression for the first Lyapunov coefficient at the Hopf 
point has been obtained in Theorem 13. 1L More concrete consequences of this 
calculation have been synthesized in Theorems 13.61 and 13.81 These results 
give sufficient conditions for the stability of the points on the Hopf hyper- 
surface and of the periodic orbit that bifurcates from the Hopf point for the 
Hexagonal governor differential equations ([6]) in two particular cases easier 
to visualize with the help of numerical plotting. See Figs. [3] and the site 

H3|. 
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